Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



Photo-ionization modelling of planetary nebulae 

II. Galactic bulge nebulae, a comparison with literature 

results 



On 



P.A.M. van Hoof^'^ and G.C. Van de Steene^'^ 



USA 



^ Kapteyn Astronomical Institute, P.O. Box 800, 9700 AV Groningen, The Netherlands 
"^University of Kentucky, Dept. of Physics and Astronomy, 111 CP Building, Lexington, KY 40506-0055, 
^European Southern Observatory, Casilla 19001, Santiago 19, Chile 

"^Research School for Astronomy and Astrophysics, Mount Stromlo Observatory, Private Bag, Weston Creek P.O., ACT 2611, Australia 



Oh! 
< 

O 



> 

(N 

O 

0\ 

o 



received, accepted 



ABSTRACT 

We have constructed photo-ionization models of five galactic bulge planetary nebu- 
lae using our automatic method which enables a fully self-consistent determination of 
the physical parameters of a planetary nebula. The models are constrained using the 
spectrum, the IRAS and radio fluxes and the angular diameter of the nebula. We also 
conducted a literature search for physical parameters determined with classical meth- 
ods for these nebulae. Comparison of the distance independent physical parameters 
with published data shows that the stellar temperatures generally are in good agree- 
ment and can be considered reliable. The literature data for the electron temperature, 
electron density and also for the abundances show a large spread, indicating that the 
use of line diagnostics is not reliable and that the accuracy of these methods needs 
to be improved. Comparison of the various abundance determinations indicates that 
the uncertainty in the electron temperature is the main source of uncertainty in the 
abundance determination. The stellar magnitudes predicted by the photo-ionization 
models are in good agreement with observed values. 

Key words: methods: data analysis - planetary nebulae: general - planetary nebu- 
lae: abundances - planetary nebulae: individual: H 1-40; M 1-20; M 2-4; M 2-23; 
M 3-15 



1 INTRODUCTION 

In van Hoof & Van de Steene (1999, Paper I) we presented 
and tested a new method to derive simultaneously and self- 
consistently all physical parameters of a planetary nebula 
from a set of observed quantities. A modified version of the 
photo-ionization code CLOUDY (Ferland 1993) is used to cal- 
culate various models, searching for a best fit of the predic- 
tions to the observables in an automated way. This method 
uses emission line ratios, the angular diameter, the radio and 
the infrared fiux to constrain the model. It also takes dust 
into account in the radiative transport. With this method we 
are able to determine the stellar temperature and luminos- 
ity, the inner, Stromgren and outer radius of the nebula, the 
density, the dust-to-gas mass ratio and the abundances. We 
investigated the accuracy of the determination of the phys- 
ical parameters by applying this method to an artificial set 
of observables. First we proved that this method can pass 
a formal convergence test. Subsequently we introduced ei- 
ther measurement errors in the observables or changed the 
model assumptions, and investigated how this affects the 



best-fit model. In this way we gained an understanding of 
the robustness of our method and hence of the reliability 
of the physical parameters. Our method was also compared 
with classical methods to determine the electron tempera- 
ture and density and nebular abundances. It was shown that 
our method suffers less from noise in the spectrum than the 
classical line diagnostics. However, this advantage may be 
lost if the model assumptions are not appropriate for the 
nebula being studied. The weakest points are currently the 
use of a blackbody approximation, the assumption that the 
inner dust radius coincides with the inner gas radius and the 
assumption of spherical symmetry. 

Distance determinations of Planetary Nebulae (PNe) 
are still very problematic. Various methods are in use, but 
the range in distances obtained is often very large and no 
method has found general acceptanc e. Rev iews of th e cur - 



rent status ca n be f ound in Pottasch (1992), Terzian ( 199. 



and Pottasch (1996). The lack of a reliable, generally appli- 
cable method to determine distances to PNe poses a prob- 
lem when using photo-ionization models. To circumvent this 
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problem we applied the method to a small sample of galactic 
bulge nebulae, which can be assumed to be all at the same 
well-known distance. Our aim is to study the accuracy of 
the determination of physical parameters by comparing our 
results with other literature values. 

A summary of the method and model assumptions is 
given in Section |2| The sample selection is presented in Sec- 
tion ^ and the modelling results in Section ^ Each PN in 
the sample is discussed individually, with special emphasis 
on the problems encountered during the modelling in Sec- 
tion ^. The resulting physical parameters are discussed by 
comparing them with results from other studies in the liter- 
ature in Section ^ Our conclusions are given in Section Q 



2 SUMMARY OF THE MODEL 

ASSUMPTIONS AND THE METHOD 

The model assumptions and the method were extensively 
described and discussed in Paper I. This section presents 
only a brief summary. 

To model the planetary nebula, we use a modified ver- 
sion of the photo-ionization code CLOUDY 84.06 (Ferland 
1993). 

The model for the PN is quite simple, and comprises 
the following assumptions: 

(1) The central star has a blackbody spectrum. 

(2) The nebula is spherically symmetric. 

(3) The density is constant inside the Stromgren radius 
of the nebula, and varies as 1/r^ outside. 

(4) Dust grains are intermixed with the gas at a constant 
dust-to-gas mass ratio; if no information on the composition 
is available they are assumed to be a mixture of graphite and 
silicates. 

(5) The filling factor, describing the small scale dumpi- 
ness of the gas, can be fixed at any value. If no information 
is available it is taken to be unity. 

(6) The distance to the nebula is fixed by an independent 
individual or statistical method. 

The above assumptions leave the following free param- 
eters: the stellar temperature, the luminosity of the central 
star, the hydrogen density in the ionized region, the inner 
radius of the nebula, the dust to gas ratio, and the abun- 
dances in the nebula. 

The outer radius of the nebula is not fixed as an input pa- 
rameter, but calculated from the long wavelength end of the 
dust emission, or, as a fail-safe, when the electron density 
drops below 0.1 cm"''. 

Adopting certain values for the input parameters, it is 
possible to calculate a model for the nebula with CLOUDY, 
predicting the continuum and line fluxes, photometric mag- 
nitudes (including the contribution of line emission) and the 
Stromgren radius. 

To compare the model predictions with the observed 
quantities, a goodness-of-flt estimator is calculated. This es- 
timator is minimized by varying all the input parameters of 
the model, using the algorithm AMOEBA (Press et al. 1986) . 

It is assumed that there exists a unique set of input pa- 
rameters, for which the resulting model predictions give the 
best fit to a given set of observables. These input parame- 



Table 1. Our sample of PNe selected from Ratag et al. (1997) 







a 


(2000) 


<5(2000) 


name 


PN G 


h 


m s 


o 1 II 


H 1-40 


359.7-02.6 


17 


55 36.0 


-30 33 33 


M 1-20 


006.1-1-08.3 


17 


28 57.5 


-19 15 53 


M 2-4 


349.8-1-04.4 


17 


01 06.2 


-34 49 39 


M 2-23 


002.2-02.7 


18 


01 42.6 


-28 25 44 


M 3-15 


006.8-1-04.1 


17 


45 31.6 


-20 58 02 



ters are then considered the best estimate for the physical 
properties of the PN. 

The full set of observed quantities necessary to derive 
the physical parameters of a PN are: 

(1) The emission line spectrum of the nebula. Usually 
this is an optical spectrum, but might also be an ultravi- 
olet and/or infrared spectrum. The line ratios make it pos- 
sible to constrain the stellar temperature, the density and 
the electron temperature in the nebula. They are also re- 
quired to determine the abundances. For elements for which 
no lines are available we assume standard abundances (AUer 
& Czyzak 1983). 

(2) Since dust is included in the model we also need in- 
formation on the mid- and far-infrared continuum. For this 
the IRAS fluxes are used. 

(3) To constrain the emission measure, either an optically 
thin radio continuum measurement (e.g. at 6 cm) is needed, 
or the absolute flux value of some hydrogen recombination 
line (usually H/3). 

(4) An accurate angular diameter of the nebula is 
needed, which we define as Od ~ 2rsti/D. Here rstr stands 
for the Stromgren radius of the nebula and D is the distance 
to the nebula. 



3 THE SAMPLE OF GALACTIC BULGE PNE 

We selected a small sample of galactic bulge nebulae from 
Ratag et al. (1997, RPDM). Galactic bulge nebulae can be 
assumed to be all at a distance of approximately 7.8 kpc 
(Feast 1987). We chose the nebulae from RPDM since they 
publish good quality spectra and also carried out their own 
photo-ionization analysis of the data which we can use for 
comparison. The radio o bserv ations for these PNe are de- 
scribed in Gathier et al. (1985). The following selection cri- 
teria were used: 



(1) The PNe should have a quality 2 or 3 IRAS 12 ^im 
fiux and quality 3 IRAS 25 /xm and 60 /xm fiuxes. 

(2) The absolute value for the radial velocity should be 
larger than 100 kms~^. 

(3) The excitation class should not be labelled peculiar. 

The resulting five PNe are presented in Table | l|. A ll nebulae 
except M 2—4 are indicated by Acker et al. (1992) as likely 
bulge PNe. In view of the large radial velocity of M 2—4, 
VhSR = —175.8 kms~^ (Gathier et al. 1983) it is unlikely to 
be a foreground object. 
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4 MODELLING RESULTS 

In Table ^ we give the input values for the observables used 
for the modelling, together with the resulting model predic- 
tions. As can be seen from this table, not all the lines present 
in the spectra are predicted by CLOUDY, most notably the 
higher Balmer lines of hydrogen and several helium lines. 
Also the element chlorine is not included in the code. The 
resulting physical parameters for the nebulae are given in 
Table ^. The hydrogen density shown in this table is the 
constant density within the Stromgren sphere. 



5 INDIVIDUAL REMARKS 

The PNe in our sample all have nearly the same medium 
excitation class. This probably is partially a result of our se- 
lection criterion that the nebulae should have been detected 
by IRAS in the 12 fim band (criterion 1). Old bulge PNe, 
having a high excitation class and cool dust, might have 
insufficient 12 fim ffux to be detected by IRAS. 

In the rest of this section each of the PNe in our sample 
will be discussed individually, with special emphasis on the 
problems encountered during the modelling. 



Table 3. The physical parameters of the galactic bulge PNe in 
our sample determined with CLOUDY. Abundances of elements for 
which only one line was observed are marked uncertain. Since we 
only model the core region of M 2—23, no values for the outer 
radius and total shell mass are entered 





H 1—40 


M 1—20 


M 2—4 


M 2—23 


M 3—15 


log{T./K) 


4.800 


4.774 


4.705 


4.782 


4.916 


log{L*/LQ) 


3.798 


3.607 


3.555 


3.639 


3.663 


log(nH/cm 


4.321 


4.124 


3.923 


4.855 


3.527 


nn/mpc 


13 


0.21 


11 


9 


33 


fstr/mpc 


24 


34 


40 


13 


98 


fout/mpc 


280: 


360: 


350: 




520: 




0.042 


0.092 


0.088 


0.015 


0.47 


M,h/MQ 


1.3: 


2.3: 


1.9: 




6.5: 


logr 


-1.70 


-3.11 


-2.50 


-2.46 


-2.60 


e(He) 


10.96 


11.02 


10.96 


11.05 


11.03 


6(N) 


7.78 


7.81 


8.13 


7.67 


7.60 


6(0) 


8.23 


8.72 


8.84 


8.67 


8.22 


e(Ne) 


7.39: 


7.82: 


8.15: 


7.75: 


7.53: 


.(S) 


6.37 


6.66 


6.90 


6.79 


6.31 


.(Ar) 


5.98 


5.99 


6.36 


6.08 


6.20: 


Te/kK 


12.7 


9.5 


8.3 


10.2 


12.0 


logf/ 


-1.40 


+2.14 


-1.16 


-1.80 


-1.58 



A colon indicates that the value is uncertain. 



5.1 H 1 40 

Two lines were omitted from the list of observables because 
of the following reasons. First the Hen A4686 line was omit- 
ted, because the flux ratio given by RPDM is quite high, 
indicative of a high stellar temperature. However, the rest of 
the observational data are not consistent with such a high 
stellar temperature. Also, this line is listed in Table 3 of 
RPDM, but is not present in their Table 1. Webster (1988, 
W88) took a spectrum of this PN, and she didn't report 
the detection of this line. She should however have detected 
a line of the strength mentioned by RPDM. Tylenda et al. 
(1994) list an upper limit of 5 for the intensity of this line. 
In view of these uncertainties we decided to omit this line. 
Since RPDM included this line in their modelling, this prob- 
ably explains the higher stellar temperature they obtain. 

The fitting of the [O ill] A4363 line was also problematic. 
The observed flux was far too low to be consistent with the 
electron temperature predicted by our model. Since the elec- 
tron temperature derived from the [N ll] line ratio is much 
higher (and more consistent with the value determined by 
our model), and also because the [Olll] A4363 line is much 
stronger in the spectrum of W88 (however not as strong as 
predicted by our model), we decided that its value was too 
uncertain and omitted it from the input. 



5.2 M 1 20 

The intensity of the Ha line seems quite high, and is not 
fitted well. The discrepancy is too large to be attributed to 
measurement errors, hence this might indicate that the spec- 
trum has not been sufficiently dereddened. There is however 
no evidence from the flts to the other lines to support this 
suspicion. 

Our model gives a very small inner radius, also resulting 
in a very high ionization parameter. This is caused by the 



high IRAS 12 fim over 25 /im flux ratio, which might indicate 
the presence of hot dust. See also the discussion in Paper I. 



5.3 M 2-4 

The spectrum is fltted quite well, but there is slight discrep- 
ancy for the [O ill] A4959 and A5007 lines. This is caused by 
the [On] A3727 doublet, which is not fitted well. The latter 
doublet usually has a larger uncertainty due to extinction 
and detector insensitivity. 



5.4 M 2 23 

This PN has the highest of all PNe in our sample. This is 
mainly caused by the weak lines, which might indicate that 
this spectrum has a lower signal-to-noise when compared to 
the other spectra. RPDM do not list error margins for their 
line fiux ratios, so we had to assume reasonable values. 

The model is not able to fit the IRAS 25 fim flux, which 
is very high compared both to the 12 fim and 60 fim flux. 
A possible explanation could be the presence of a 30 fim 
dust feature in the spectrum (Hoare 1990). This would imply 
that the nebula is carbon-rich, since this feature has only 
been observed in carbon-rich nebulae. The central star has 
spectral type Of (AUer & Keyes 1987, AK87). 

The large difference between the optical diameter of 
8.5" (Acker et al. 1992) and the radio diameter of 0.72" 
(Gathier et al. 1983) suggests that this nebula might be a 
core-halo nebula. All other evidence gathered in this paper 
also is consistent with this assumption and we will adopt it 
throughout the paper. Since we used the radio diameter for 
the modelling, our model is only valid for the core region. 
The fact that our model is density bounded and gives a low 
ionized mass is consistent with the fact that we are only 
modelling the core region. 
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Table 2. Comparison of the observed quantities (mostly taken from Ratag et al. 1997) and the model fit for our sample of PNe. The 
strength of the emission lines is given relative to H/3 = 100. The measured line fluxes have been dereddened. The entries [O ll] A3727, 
A7325 and [So] A4071 all stand for the entire multiplet. All observables for which entries in both columns obs. and model are present, 
have been weighted in the goodness-of-fit estimator, except where indicated. 



ion 


A 


H 1 


-40 


M 1 


-20 


M 2-4 


M 2 


-23 


M 3- 


-15 




A 


obs. 


model 


obs. 


model 


obs. 


model 


obs. 


model 


obs. 


model 


[On] 


3727 


32.4: 


38.5 


55.5 


83.1 


97.1 


125.2 


14.2 


20.4 


48.6 


101.5 


H12 


3750 














3.4 








Hll,Om 


3771 














4.2 








HIO 


3798 






5.5 




4.0 




5.1 








H9 


3835 






6.7 




7.1 




6.5 








[Nciii] 


3869 


79.1 


79.4 


69.1 


67.6 


64.7 


67.6 


82.4 


80.0 


89.9 


90.1 


H8,Hei 


3889 


12.9 




18.1 




17.4 




13.9 




16.6 




[Neni],He 


3969 


19.9 




25.5 




33.6 




38.3 




25.4 




He I 


4026 






2.7 




2.7 




1.82 








[Sii] 


4071 


3.3 


2.5 




2.3 


3.6 


3.0 


2.5 


3.8 




1.8 


H5,Nin 


4102 


25.8 


29.4 


25.6 


30.3 


24.7 


30.5 


24.0 


30.2 


26.1 


29.5 


Cii 


4267 






















H7 


4340 


47.1 


50.2 


49.1 


51.1 


45.9 


51.2 


48.2 


50.9 


49.0 


50.2 


[Oiii] 


4363 


4.6* 


12.0 


7.5 


5.4 


2.9 


2.1 


13.9 


9.9 


3.3: 


7.5 


He I 


4472 


6.5 


4.4 


5.8 


5.2 




4.5 


5.1 


5.5 


5.0 


5.1 


Niii 


4641 


















20.3 




Hen 


4686 


17.1?t 


1.9 




0.4 




0.14 




0.6 


4.1 


3.5 


[Ar iv] ,He I 


4712 






1.33 




0.91 




1.09 




1.94 




[Ar iv] 


4740 


4.45 


4.70 


0.71 


1.16 


0.33: 


0.69 


0.80 


1.42 




6.67 


H/3 


4861 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


100. 


He I 


4922 


1.36: 








1.42 




0.87 








[Oiii] 


4959 


307. 


276. 


336.t 


315. 


272.t 


213. 


304. 


350. 


328. t 


219. 


[Oiii] 


5007 


915. 


827. 


1009.t 


946. 


815.t 


640. 


1006. 


1051. 


983. t 


658. 


[Ni] 


5201 




0.20 




0.07 




0.22 




0.00 




0.59 


[CI hi] 


5517 














0.21 








[CI hi] 


5538 






0.47 








0.25 








[Nil] 


5755 


1.92: 


2.16 


0.92 


1.05 


1.39 


1.52 


1.19 


1.06 


1.0: 


1.4 


He I 


5876 


15.5 


15.9 


15.7 


16.0 


13.8 


13.8 


17.2 


17.8 


16.3 


16.3 


[Oi] 


6300 


2.8 


2.7 


4.7 


4.1 


2.9 


3.7 


4.2 


3.3 


2.8 


3.6 


fSiiil 


6312 


1.70 


1.48 


0.79: 


1.12 


1.59 


1.33 


2.4 


2.3 


1.24 


1.00 


[Oi] 


6364 


0.88 


0.88 


1.46 


1.35 


1.00 


1.23 


1.42 


1.09 


0.79 


1.19 




6548 




19.9 




14.6 




32.1 




6.4 




18.8 


Ha 


6563 


280. 


278. 


303. 


269. 


275. 


268. 


283. 


269. 


305. 


278. 


[Nil] 


6584 


61.4 


59.7 


45.4 


43.8 


85.6 


96.2 


18.9 


19.2 


57.6 


56.3 


He I 


6678 


3.8 




4.0 




3.2 




4.7 




4.6 




[Sii] 


6716 


0.99 


0.96 


1.17 


1.45 


2.72 


3.16 


0.78 


0.47 


2.65 


3.38 


[Sii] 


6731 


1.71 


1.99 


2.32 


2.91 


5.0 


6.0 


1.55 


1.05 


5.4 


5.2 


He I 


7065 


7.9 


11.4 


10.5 


9.7 


5.6 


7.1 


14.4 


12.0 


7.5 


9.2 


[Ar III] 


7136 


13.6 


13.3 


9.1 


7.5 


15.2 


13.0 


14.0 


11.2 


19.2 


19.3 


He I 


7281 






0.83 




0.3: 




0.94 




0.59 






7325 


9.2 


14.2 


14.9 


17.1 


8.3 


14.0 


19.3 


21.4 


6.6 


8.8 


obs. 


unit 






















F^(12 fim) 


Jy 


2.38 


2.38 


1.13 


1.00 


0.56 


0.53 


1.93 


2.10 


< 0.53§ 


0.19 


F„ {25 fim) 


Jy 


18.45 


19.11 


3.94 


4.44 


5.00 


5.83 


9.31 


6.54 


5.66 


6.02 


F^(60 Aim) 


Jy 


11.91 


11.42 


2.38 


2.30 


5.77 


5.18 


1.64 


1.64 


8.02 


7.77 


F^(100 Atm) 


Jy 


< 73.48 


3.22 


< 4.59 


0.66 


< 12.59 


1.73 


< 126.70 


0.24 


< 10.39 


2.72 


F^(6 cm) 


mJy 


31. 


31.0 


47. 


47.7 


32. 


32.2 


41. 


41.5 


65. 


65.4 


ed 


arcsec 


1.26 


1.27 


1.98 


1.81 


2.16 


2.13 


0.72 


0.67 


5.4 


5.19 








0.63 




1.68 




2.28 




4.17 




2.75 



: A colon indicates that the value is uncertain. 

t The sum of the intensities of the doublet was split using the ratio 3:1. 

■t This line was not weighted in the goodness-of-fit estimator , see also Section ^ 

§ This flux is not listed as an upper limit in the IRAS Point Source Catalog, see also Section bl 
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5.5 M 3 15 

There is a suggestion of a systematic trend when comparing 
the observed and the modeUed Une flux as a function of 
wavelength. Also the observed intensity of the Ha line seems 
quite high. This might indicate that the spectrum has not 
been sufficiently dereddened. 

This PN has a [WC]-type central star (AK87). The cen- 
tral star temperature, the excitation class 5.5 (taken from 
RPDM) and the low IRAS 12 fim to 25 fim flux ratio all are 
consistent with an early spectral type: [WC3-4] (cf. Kaler 
1989, Mendez & Niemela 1982, and Zijlstra et al. 1994, re- 
spectively). 

The IRAS 12 /im flux is not listed as an upper limit 
in the Point Source Catalogue. However, when we used this 
value for the modelling, the resulting model was unrealistic. 
We therefore assume that the 12 nm flux suffers from con- 
fusion and took the quoted value to be an upper limit. See 
also the discussion in Paper I. 



6 DISCUSSION 

In this section the modelling results are discussed by com- 
paring them with results from other studies in the literature. 
Since distance dependent parameters are usually not given 
by other authors, we will restrict ourselves to the distance 
independent parameters of PNe. 



6.1 Stellar temperatures 

In Table ^ we present a comparison with the stellar temper- 
atures given in the literature. These temperatures were de- 
rived using the Zanstra method and photo-ionization mod- 
elling. Results using the energy balance or Stoy method are 
not listed since we consider this method, or at least the 
data for the nebulae being studied here, to be unreliable 
(Pottasch, private communication). One can see that the 
derived values agree quite well with only a few outliers. 

The temperatures determined by our method agree well 
with the hydrogen Zanstra temperatures, with the single 
excepti on of t he temperature for M 2—23 given by Tylenda 
et al. (1991a). Since the other three determinations using 
the Zanstra method agr ee quit e well, we assume the value 
given by Tylenda et al. (1991a) to be erroneous. 

To derive stellar temperatures for photo-ionization 
modelling, sometimes certain line-ratios are used as temper- 
ature indicators (e.g. Hen A4686 over H/3). Especially for 
cooler central stars, where few temperature sensitive lines 
are available, this makes the determination dependent on 
one or two lines. Nevertheless, the results from the photo- 
ionization models usually are in good agreement. Exceptions 
are the temperature for H 1—40 derived by RPDM, and the 
temperatures for M 3—15. The deviating value for H 1—40 
given by RPDM can probably be attributed to the Hen 
A4686 line, which they used as a temperature indicator. We 
refer to the discussion in Section 5.1, For M 3—15 we find a 
higher stellar temperature than other authors. The largest 
discrepancy is with the value from AK87. This can prob- 
ably be attributed to the fact that AK87 did not report 
a detection of the Hell A4686 line in their spectrum (al- 
though a detection of roughly the same strength as RPDM 



Table 4. Comparison of the stellar temperatures for the PNe in 
our sample. The temperatures are given in kilokelvin. The abbre- 
viations for the methods have the following meaning: H I - hydro- 
gen Zanstra method. He II - helium Zanstra method, AM - photo- 
ionization modelling using model atmospheres, BB - photo- 
ionization modelling using blackbody approximation. 



1-40 


M 1-20 


M 2-4 


M 2-23 


M 3-15 


ref. 


meth. 






55. 


64. 




3 


Hi 




53. 




56. 




4 


Hi 






51.5 


65.0 




6 


Hi 




65. 




85. 




8 


Hi 




49.5 


49.9 






5 


Hi 




59.9 








5 


He n 








50. 


62.5 


1 


AM 


80.0? 


50.0 


50.0 


57.5 


72.5 


7 


AM 


64. 










2 


BB 


63.1 


59.4 


50.7 


60.5 


82.4 


9 


BB 



References — 1. Aller fc Keyes (1987) us ing m odel atmospheres 
by Husfeld et al. (^984 ) 2. Dopita et al. f|l99o| ) 3. Gleizes, Acker 
& Stenholm (|l989|) 4. Kaler & Jacoby (|l"99l|) 5. Mal'kov (|l997|) 



6. Pottasch & Acker ( |l989[ ) 7. Ratag et al. ( |l997| ) using model 
atmos pheres by Clegg & Midd lemass (1987) and Husfeld et al. 



1984) 8. Tylenda et al. ( L991a ) 9. This work 



was reported in Aller & Keyes 1985). Since M 3—15 has 
a [WC]-type central star, part of the Hen A4686 fiux may 
originate from the central star. Unfortunately, no detection 
of the [Ariv] A4740 line has ever been reported, so that no 
alternative temperature sensitive line is available. In view 
of this, the central star temperature for M 3—15 should be 
viewed with some caution. 

We conclude that the temperature determination for 
the central stars in this sample is fairly reliable, although the 
situation for M 3—15 is not completely clear. This confirms 
our results from Paper I in which we found the temperature 
determination to be robust. 



6.2 Electron temperatures 

In Table |^ the electron temperatures derived by different 
authors are compared. The electron temperature determined 
by CLOUDY is a weighted mean of the temperature in the 
nebula: Te = J n^TcdV/ j n^dV. The observational material 
shows a large spread in most cases, even when the same 
method is used. This indicates that the electron temperature 
determination, at least in those cases where diagnostic lines 
have been used, is not very reliable. This is in agreement 
with our results in Paper I. Note the large difference between 
the [N ll] and [O ill] temperatures in the case of M 2—23. 
This difference is not caused by measurement error. For this 
particular object, the temperature derived from the [Nil] 
lines has no physical meaning (Liu, private communication) . 

The electron temperatures derived from our method are 
in most cases just outside the range of values found with 
line diagnostics; three times at the low end and twice at 
the high end. The results from Paper I indicate that the 
electron temperature determination with our method should 
be robust. It is not apparent to us why the average values 
of the electron temperature derived from line diagnostics do 
not coincide with our results. This might indicate a problem, 
although the fact that we find both higher and lower results 
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Table 5. Various determinations of the electron temperature 
for the nebulae in our sample. The temperatures are given in 
kilokelvin. The abbreviations for the methods have the following 
meaning: ave. - average of [N ll] and [O III] , model - average model 
prediction (see text). 



Table 6. Various determinations of the electron density for the 
nebulae in our sample. The densities are given in 10'' cm~^. The 
abbreviations for the methods have the following meaning: radio - 
density determined from the radio flux, model — average model 
prediction (see text). 



H 1-40 


M 1-20 


M 2-4 


M 2-23 M 3-15 


ref. 


meth. 








18.0 


1 


[N II] 




17.1 






2 


[N II] 




12.5 


12.5 




4 


[Nil] 






10.2 




7 


[Nil] 




10.4 






8 


[Nil] 


13.1 


10.2 


9.7 


19.2 9.4: 


9 


[N 11 


10.1 








11 


[Nil] 








11.0 


1 


[0 in] 




10.8 






2 


[0 III] 




13.1 


11.6 




4 


[O III] 






8.7 




5 


[O III] 








12.9 


6 


[O III] 




9.9 






8 


[O III] 


9.3 


10.4 


8.5 


13.0 8.4: 


9 


[O III] 


9.7 








11 


[O III] 








12.6 11.2: 


3 


ave. 






11.1 




10 


ave. 


12.7 


9.5 


8.3 


10.2 12.0 


12 


model 



: A colon indicates that the va 
References — 1 . Ac ker et al 
AUer 



uc is uncertain. 
( |1989| 2. A cker et al. ( 1991 ) 
Kevcs (|1987|) 4. Costa et al. |l99(^) 5. Cuisin icr. Acker , 



Kop pen (|1996| ) 6. Kaler ( |l97q) 7 . Kaler et al. ( |l993| ) 8. Kal cr et 
al. C 1996 ) 9. R atag et al. ([[997[) 10. Tylenda et al. (|l9911:|) 11. 



Webster (1988) 12. This work 



is not indicative of a systematic effect. Nevertheless, this 
issue should be investigated further in future research, using 
a larger sample. 



6.3 Electron densities 

In Table ^ the electron densities derived by different au- 
thors are compared. The electron density determined by 
CLOUDY is a weighted mean of the density in the nebula: 
= J n^dV/ J n^dV. There are enormous differences between 
the various determinations in the literature, even when the 
same method has been applied. This indicates that the de- 
termination of densities with line diagnostics is unreliable, 
which confirms our results in Paper I. Also note the enor- 
mous differences between the [S ll] , [CI il l] and [Ar iv] densi- 
ties for M 1-20 derived by Kaler et al. (|l996|). 



Our values differ substantially from the values given by 
RPDM, although they are based on the same observational 
data. This is because we use a completely different method 
to determine the density. For three out of five nebulae we 
find results which are within the range of values found with 
other methods. For M 2—4 we find a value which is a bit 
larger. The results in Paper I indicate that our determina- 
tion of the density is somewhat susceptible to measurement 
errors and errors in the model assumptions. This might pro- 
vide an explanation for the discrepancy. The fact that we 
model only the core region of M 2—23 provides an expla- 
nation fo r the very high density we find for this nebula. 
Webster ( |l976| ) and Bofii & Stanghelhni (1994, using the 
same spectrum) also find a high value using the [Ar iv] line 



H 1-40 


M 1-20 


M 2-4 


M 2-23 M 3-15 


ref. 


meth. 








13.6 


1 


[S II] 


10.9 


4.7 




24.3 


2 


[Sii] 








3.0 2.5 


3 


[Sii] 




17.8 


4.5 




5 


[S II] 






7.0 




6 


[S II] 


15.0 








7 


[S II] 






4.2 




9 


[S II] 




85. 






10 


[S ii] 


4.4 


9.2 


5.6 


11.5 10.6 


11 


[S II] 








4.2 


13 


[S II] 






3.6 




14 


[S II] 


35.1 








16 


[Sii] 






5.7 




9 


[CI III] 




7.8 






10 


[CI III] 








79. 


4 


[Ar iv] 




1.0 






10 


[Ar iv] 








63. 


15 


[Ar iv] 


13.5 








7 


radio 








20. 


8 


radio 




10. 






12 


radio 


22.7 


14.6 


9.1 


79.3 3.7 


17 


model 



References — 1 Acker et al. ( |l989| ) 2. Ack er et a l. ( |l99l| ) 3. AUer 
fc Kc vcs (|l987|) 4. BofH fc St anghellini (|l994|) 5. Co sta et al. 



1996| ) 6. Cuisinicr et al 
( |1979| ) 9. K aler et al. ( [l993| 



1996| ) 7. Dopita ct al, 
Kalcr ct al. 



1990) 



10. 



Kaler 
11. Ratag 

Kalcr (|1989|) 13. St anghe llini & Kaler 



et al. ( |1997| ) 12. Shaw 

( 198£ ) 14. Tylenda et al. (|l991b|) 15. Webster ( |l97l 
(1985) 17. This work 



16. Webster 



ratio. The [Ariv] lines are expected to be formed predomi- 
nantly in the core region and hence this would confirm our 
results. On the other hand, the excitation in the core is too 
high to form large amounts of S""" . Ifence the [S ll] lines can 
be expected to originate predominantly from the halo and 
should therefore indicate lower densities. This of course also 
depends on the exact position of the slit over the nebula. 
All of this might be an explanation for the extremely large 
spread of values found for this nebula. 

The quality of the data in Table ^ makes a comparison 
with our results meaningless. However, the data are at least 
consistent with the assumption that our results are more 
accurate than the results from line diagnostics. 

6.4 Nebular abundances 

In Table ^ we give a comparison of the abundances we de- 
termined with other literature values. We did not include 



the nitrogen abundance for M 3—15 from Henry (199C). Af- 
ter a discussion with Dr. Henry it was established that this 
abundance was flawed by an error in the analysis (as is also 
the case for the nitrogen abundances of M 4—3 and H 1 — 23 
listed in the same paper; all other results are not affected). 
We also did not include the abundances for M 2—23 listed 



in Koppen, Acker & Stenholm (1991). It was established 
that this analysis was flawed by an error as well, and Dr. 
Koppen kindly provided us with a re-analysis of his data. 
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Table 7. Comparison of the abundance determinations of the PNe in our sample. 







H 1- 


-40 








M 1- 


-20 








M 2 


-4 




ref: 


4 


9 


12 


13 




2 


7 


9 


13 




2 


3 


9 


13 


e(He) 


11.03 


11.06 


11.04 


10.96 




11.07 


10.94 


11.02 


11.02 




11.11 


10.96: 


10.99 


10.96 


e{N) 


7.72 


8.08 




7.78 




7.39 


7.80 


7.75 


7.81 




7.65 


8.17 


8.09 


8.13 




8.52 


8.70 


8.53 


8.23 




8.30 


8.65 


8.62 


8.72 




8.30 


8.77 


8.80 


8.84 




7.89 


7.69 




7.39: 








7.79 


7.82: 








7.90 


8.15: 




6.88 


6.77 




6.37 




6.43: 




6 52 


6 66 




6.64 


7.03 


6.96 


6.90 




6.6 


6.43 




5.98 








6.05 


5.99 




6.31 


6.36 


6.25 


6.36 










M 2- 


-23 












M 3- 


-15 






ref: 


1 


5 


6 


8 


9 


10 


11 


13 




1 


5 


9 


13 




e(He) 


11.00 


10.93 




10.88: 


10.98 


10.92 


10.96 


11.05 




11.03 


11.01 


11.03 


11.03 




e(N) 


7.68 


8.20 


7.70 


7.55 


7.40 


8.13 




7.67 




8.08 




8.14 


7.60 




6(0) 


8.40 


8.42 


8.18 


8.34 


8.22 


8.47 


8.11 


8.67 




8.41 


8.51 


8.74 


8.22 




e(Ne) 


7.60 


7.62 




6.46: 


7.15 


7.65 




7.75: 




7.48 


7.62 


7.86 


7.53: 




6(S) 


6.6 






6.29 


6.30 






6.79 




6.7: 




6.86 


6.31 




e{Ar) 


5.75 






5.86 


5.81 






6.08 




6.5 




6.53 


6.20: 





: A colon indicates that the val ue is uncertain. 

Refe rences — 1. AUer fc K cycs (1987) 2. Costa et al. (1996) 3. Cuisinier et al, 
( | l99 6[) 



( ) 7. Kaler et al 
(1976) 12. Webster 



( 1996 ) 4. Dopita et al 
Koppen (private communication) 9. Ratag et al. (1997) 10. Walton, Bar 



(|l99q) 5. Hen ry (|l990[ ) 6. Kaler 
z Clegg (|l9935ll. Webster 



13. This work 



The higher nitrogen abundance given by Walton et al. ( 199; ) 
for M 2—23 might be a result of the inclusion of lUE data 
in their analysis. They systematically find higher nitrogen 
abundances for bulge PNe than other authors. 

One can see that large differences can be found between 
the various abundance determinations in the literature. If we 
exclude our own results, we find the following statistics. For 
elements heavier than helium, we find a difference between 
the lowest and highest abundance determination larger than 
or equal to 0.3 dex in 12 out of 22 cases, and larger than 
or equal to 0.5 dex in 4 out of 22 cases. For the helium 
abundances we find a spread larger than 0.1 dex in 2 out 
of 5 cases. Especially the abundances for M 2—23 show a 
large spread and should be considered uncertain. From this 
we draw the conclusion that, at least for the sample studied 
here, abundance determinations can not be considered very 
accurate. Uncertainties exceeding 0.2 dex to 0.3 dex are not 
uncommon. 

When one compares the abundances for the individual 
PNe with the values from RPDM, one can see that for the 
two objects where the electron temperature is in good agree- 
ment (M 1—20 and M 2—4), the abundances also agree very 
well. For the other objects the abundance determinations 
differ. We attribute this to the difference in the determi- 
nation of the electron temperature. When we compare our 
abundance determinations with the other values found in the 
literature, we see that our results often are slightly outside 
the range of values found by other authors. This behaviour is 
well correlated: either all outliers are at the low end or at the 
high end. This behaviour is also well correlated with our elec- 
tron temperature determination. When our electron temper- 
ature determination is at the low end, our abundances are at 
the high end, and the reverse is also the case (see also Sec- 
tion 3.2). This indicates that the main source of uncertainty 



in the abundance determination is th e ele ctron temperature. 
Hence the discussion given in Section 3.2 applies here as well. 



6.5 Stellar broadband fluxes 

Since CLOUDY calculates the attenuation of the stellar con- 
tinuum separately from the transport of the diffuse nebular 
continuum, we are able to predict broadband photometric 
fluxes for the central star as they would appear through 
the nebula. In this way we could calculate a prediction for 
the Johnson B and V magnitudes. However, observed stellar 
magnitudes will be reddened due to interstellar extinction 
as well, and we have to take this into account in our predic- 
tions. To calculate the total extinction towards the nebula, 
we averaged all the measurements we could find in the lit- 
erature. Since the continuum fluxes predicted by CLOUDY 
already take the internal extinction into account, we only 
have to correct the stellar magnitudes for the external ex- 
tinction. Hence we used the internal extinction from our 
model, and subtracted it from the total extinction. Then 
we used this value for the external extinction to predict the 
reddened B and V magnitudes of the central star. Where 
necessary, we appl ied the interstellar reddening law given 
by Pottasch ( 1984 ). A comparison of the calculate d val ues 
with the literature values taken from Acker et al. (1992) is 
given in Table ^ 

The predicted magnitudes are slightly fainter than ob- 
served, but still in remarkable good agreement, considering 
the fact that we use a blackbody approximation to deter- 
mine these values. Given the fact that a blackbody of a 
given temperature has more ionizing photons than a realis- 
tic spectrum with the same effective temperature, one can 
expect that in the best-flt model the total luminosity will 
be underestimated to compensate for this effect. However, 
we flnd that this effect is only very modest and this can be 
understood from the fact that we include the dust emission 
in the modelling. Grains can be heated very efficiently by 
Balmer continuum photons, as well as by Lyman continuum 
photons. Therefore, the IRAS fluxes give a good constraint 
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Table 8. For the PNe in our sample we give in column 2 and 3 the Johnson B and V magnitudes resp. predicted by our model, in 
column 4 the internal extinction in the Johnson V band derived from our model, in column 5 the average total extinction derived from 
the Balmer decrement and the radio flux, in columns 6 and 7 the predi cted reddened values for the Johnson B and V magnitudes and 
in column 8 and 9 the measured magnitudes given in Acker et al. (11993). 



name 


^mod 

mag 


^mod 

mag 


mag 


^tot 

mag 


^pred 

mag 


^prcd 

mag 


B 
mag 


V 
mag 


H 1-40 
M 1-20 
M 2-4 
M 2-23 
M 3-15 


15.70 
14.62 
14.41 
14.63 
15.52 


15.68 
14.90 
14.64 
14.90 
15.80 


1.11 
0.06 
0.14 
0.08 
0.10 


5.05±0.26 
2.61±0.10 
2.78±0.15 
1.81±0.22 
4.69±0.25 


20.86±0.34 
17.96±0.13 
17.87±0.19 
16.90±0.28 
21.53±0.33 


19.62±0.26 
17.45±0.10 
17.28±0.15 
16.63±0.22 
20.39±0.25 


17.7 
17.6 
16.7 


17.1 
17.0 



on the Balmer continuum flux. This counteracts the previ- 
ously mentioned underestimation of the total luminosity and 
explains the remarkable accuracy of our stellar broadband 
fluxes. 



6.6 Distances 

In our model assumptions we assume the distance to be a 
flxed number. However, our method can easily be changed 
in such a way that the distance would be a free parameter. 
When this is done, the best-fit model would also give an 
estimate for the distance. We have investigated the possibil- 
ity to determine the distance this way (van Hoof & Van de 
Steene 1996). We found that, though possible in principle, 
the spread in the resulting distance determinations is large. 
The distance determination is vulnerable to various observa- 
tional errors, but especially to the error in the determination 
of the angular diameter. Since the angular diameter is no- 
toriously hard to measure, this sensitivity makes the results 
very uncertain. When we determined the distances to the 
bulge nebulae in our sample with this method, we found the 
spread in the values to be larger than what is obtained from 
a statistical method (Van de Steene & Zijlstra 1995). 

Closer investigation reveals that this method of deter- 
mining distances is in essence id entica l to the method de- 
scribed in Phillips & Pottasch (1984). They already con- 
cluded that this method is unreliable. The use of a wrong 



value for the distance not only influences the distance de- 
pendent parameters but also some distance independent pa- 
rameters, as was already discussed in Paper I. We therefore 
do not recommend this method, and advise the use of sep- 
arately determined distances. 



7 CONCLUSIONS 

We applied our method which enables a fully self-consistent 
determination of the physical parameters of a PN, using the 
spectrum, the IRAS and radio fluxes and the angular diam- 
eter of the nebula, to a sample of five galactic bulge PNe. 
Comparison of the distance independent physical parame- 
ters with published data shows that the stellar tempera- 
tures generally are in good agreement and can be considered 
reliable. The literature data for the electron temperature, 
electron density and also for the abundances show a large 
spread, indicating that the use of line diagnostics is not reli- 
able. Comparison of the various abundance determinations 



indicates that the uncertainty in the electron temperature 
is the main source of uncertainty in the abundance deter- 
mination. The large spread in the literature data makes a 
comparison with our results meaningless. The stellar magni- 
tudes predicted by the photo-ionization models are in good 
agreement with observed values. 
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